Single hole dynamics in the one dimensional t-J model 
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We present a new finite-temperature quantum Monte Carlo algorithm to 
compute imaginary-time Green functions for a single hole in the t-J model 
on non-frustrated lattices. Spectral functions are then obtained with the 
Maximum Entropy method. Simulations of the one-dimensional case show 
that a simple charge-spin separation Ansatz is able to describe the overall 
features of the spectral function over the whole energy range for values of J/t 
from 1/3 to 4. This includes the bandwidth ~ 4t + J and the compact 
support of the spectral function. The quasiparticle weight Zk is computed on 
lattices up to L = 96 sites, and scales as Z^, oc L~^/^. 

Understanding single hole dynamics in quantum antiferromagnets is a decisive step to- 
wards a comprehensive description of elementary excitations in strongly correlated systems. 
Experimental realizations are found in compounds such as SrCu02 |Q, ^^^0.96^205 for 
chains, Sri4^Cu2404i for ladders and Sr2Cu02Cl2 @] for planes. In particular chain com- 
pounds attract at present an increasing amount of interest in order to elucidate, whether 
signals of charge-spin separation as predicted from Luttinger-liquid theory can be observed 
experimentally. On the other hand, theoretical treatments based on Bethe-Ansatz (BA) re- 
sults lead recently to a complete description of the spectral function of the Hubbard model 
at f/ = 00 1^ and the low energy sector in the nearest-neighbour (NN) t-J model, where 
explicit results are obtained at the supersymmetric (SuSy) point 0. Further exact results - 
apart from exact diagonalizations which suffer from strong finite-size effects - are available 
only for the inverse-square exchange (ISE) t-J model at the SuSy point. In order to be 
able to compare with experiments, it is crucial to extend such studies to realistic values of 
the parameters and possibly beyond the asymptotic low energy limit. 

In this letter, we present a simple finite-temperature quantum Monte Carlo (QMC) al- 
gorithm capable of dealing with this issue for the NN t-J model. For single-hole excitations 
and in the absence of frustration, the method is free of the notorious sign problem, and ap- 
plicable to chains, n-leg ladders and planes. Here, we concentrate on chains. Our simulations 
lead to the conclusion that the overall features of the spectral functions are well described 
by a charge-spin separation Ansatz (CSSA) based on a mean-field slave-boson picture 0, 
where the hole spectral function is given by the convolution of the spectral functions of free 
holons and spinous. The agreement with the simulations is obtained over all energy scales 
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and values of J/t ranging from 1/3 to 4. At the SuSy point a more detailed understanding of 
the spectrum is achieved by supplementing the simple model with BA results. A finite-size 
scaling on chains up to L = 96 sites shows that the quasi-particle weight Zk vanishes as 
1 / V^, a result which was beyond numerical capabilities up to now. 
Our starting point is the NN t-J model, 

Ht-j = -t ^Ucj,a + 'J (^i- 4- - J^i^j) ■ (1) 
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transformation this model is cast into the form 10 



H,.j =+tY: p^,fif, + 1 E ^vip^, - 1), (2) 
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where Pij = (1 + CTj ■ c?j)/2, Aij = (1 — — Uj) and rii = f-f^. In this mapping, one 
uses the following identities for the standard creation (c| g.) and annihilation (Cj ^) operators 
4^ = ji^+fi-ji-fl , cj^ = ai-{fi + f}) , where 7i,± = {l±(Ti^^)/2 and ai,± = {ai^^±iai,y) /2. 
The spinless fermion operators fulfill the canonical anticommutation relations {//, fj} = 6ij, 
and ai^a , a = x,y, or z are the Pauli matrices. The constraint to avoid doubly occupied 
states transforms to the conserved and holonomic constraint I]j7i,-///j = 0. 
The Green function in the spin up sector may be written as 

G^{z-j,t) = (Tc^^^irycl) = {Tfhr)f,) (3) 

where T corresponds to the time ordering operator. Inserting complete sets of spin states 
the quantity above transforms as 
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= J2Pi^)GihJ,T,a) + OiAr') (4) 

Here mAr = (3, nAr = r. Art ^ 1 and exp(— AriJ(o"i, (J2)) is the evolution operator 
for the holes, given the spin configuration ((Ti,o"2). In the case of single hole dynamics 
\v) is the vacuum state for holes, and -P(c?) is the probability distribution of a Heisenberg 
antiferromagnet for the configuration a, where a is a vector containing all intermediate states 
(cTi, . . . cr„, . . . (Tm, CTi). The sum over spins is performed in a very efficient way by using a 
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world- line cluster-algorithm |TT|. As the evolution operator for the holes is a bilinear form 
in the fermion operators, G{x, r, a) can be calculated exactly. G{x, r, a) contains a sum over 
all possible fermion paths between {i, 0) and (j, r), where i — j = x. This stands in contrast 
to the worm approach |T^, where fermion paths are sampled stochastically. The numerical 



effort to calculate G{x, r, a)"^ x,t scales as Lr. Spectral properties are obtained by inverting 
the spectral theorem 

oo 

G{k,T)= [ dujA{k,u)—^p^p—- (5) 
J TTi 1 + expi — pa;)) 

— oo 

with the Maximum Entropy method (MEM) Since -P(ct) is the probability distribution 
for the quantum antiferromagnet, the algorithm does not suffer from sign problems on 
bipartite lattices and next neighbour interactions in any dimension. However, when the 
spin and charge dynamics evolve according to very different time scales (J^ 0.2t), G{x, r, a) 
shows an increasing variance. Best results are obtained at the SuSy point and an appreciable 
range of J/t may be considered as shown below. 

We now concentrate on the one-dimensional t-J model. The simulations were performed 
at temperatures T < min(J, t)/15, such that no appreciable changes with a further decrease 
in temperature can be seen: the results correspond to the zero temperature limit, a limit 
which is in general difficult to reach in other finite-temperature fermionic algorithms. We 
compare our results with the predictions of the CSSA, where free holons and spinous are 
described by PJI^ 



H = -jT. hjh^ - ^ E 4,.s,,- (6) 

<i,j> <i,j> 

Here the electron operator c^ ^ is given by the product of a holon (hi) and a spinon (sj^o-) 
operator, c^^ = Si^^hl, the holon being a boson and the spinon a spin-1/2 fermion. As a 
consequence of the above Ansatz , the dispersion relations of the free holons and spinous 
are given by eh = —th cos qh and = — J^cosg^ respectively, whereas the energy of the hole 
is E{k) = eh — €s and by momentum conservation k = qh — qs- We take th and Js as two 
free parameters in contrast to a mean-field approximation, where they have to be calculated 
self-consistently. The spectral function is then given by a convolution of the spinon and 
holon Green functions. The lowest attainable energy {—th) and highest one {th + Js) define 
the bandwidth of the hole, 2th + Js- Since the full band- width obtained by considering the 
compact support of the spectral function at J = is known to be exactly 4t 0, we take 
th = 2t. In order to determine J^, we consider the overall bandwidth, as obtained from the 
simulation. As can be seen in Fig. 1, for all values of J, 
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FIG. 1. Density of states N{uj) for different values of J jt. The vertical line indicates 4i + J. 

the width of the density of states N{uj) scales approximately as 4t + J in the parameter 
range considered, leading to Jg = J- 

Beyond predicting bandwidths, the CSSA describes accurately the support of the spec- 
tral function in the case J = 0, when compared with exact results PJ^. If furthermore 
phase string effects |0 are taken into account, the singularities of A{k,u) related to holons 
and spinous can be reproduced. For finite J, the minimal (maximal) possible energy of 
a hole in CSSA is given by E{k) = —F^ {E{k) = F^) for k < {k > k^), where 
Fk = \J .P + 4t^ — 4tJ cos {k) contains both holon and spinon contributions, and k^ is de- 
termined by cos(A;o) = J/{2t). The remaining parts of the compact support are given by 
E{k) = =i=2tsin(A;) for k > ko (lower edge) and k < ko (upper edge) respectively. Such 
dispersions correspond to holons with momentum k + qs, and a spinon with = ^71/2 
As J — i> 2t, /cq — ^ and the lower edge of the compact support is entirely determined by the 
dispersion of the holon. 
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FIG. 2. Spectral function A{k,uj) for J = OAt (a), J = 1.2t (b) and J = 2t (c). Here the wave 
vector k G [0, vr] is given on the y-axis. For clarity, the data is rescaled by the number given at the 
right hand side of the plot. Further details are discussed in the text. 



We now compare the above predictions with our QMC data. Figure 2 shows A{k^uj) for 
J/t = 0.4 (a), 1.2 (b) and 2 (c). In all cases the compact support is reproduced very well 
by the CSSA. The Ansatz also predicts singularities at the lower (upper) edge for k < ko 
{k > fco), and when phase strings are considered along the edges and the holon lines 
(±2tsin(A;)) for all momenta. The singularities along the lower holon line are also supported 
by a recent low energy theory |^ . For all parameter values we observe dominant weight along 
the above mentioned lines. For J/t = 0.4, we have checked that the results are consistent 
within the uncertainties of MEM with a peak along the edges and a further peak along 
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the holon lines, signaled by a broad structure between the edges and the holon lines (Fig. 
2. a). We observed such a behaviour for 0.33 < J/t < 0.6. For J/t > 1.2 (Fig. 2.b and 
2.c), the structure at the lower edge narrows considerably and the data are not any more 
consistent with an additional structure along the lower holon line for k < k^, but only with a 
singularity for k > k^. At J/t = 2 the exact holon and spinon dispersions can be obtained by 
BA [^. Figure 2.c shows the comparison with the CSSA, where on the one side the original 
dispersions are used (full line) and on the other side, with the dispersions as given by BA 
(crosses). Whereas the BA holon dispersion reproduces very well the lower edge, showing 
that as anticipated by the CSSA, at the SuSy point that edge is completely determined by 
the holon dispersion, the full bandwidth is better described with the original dispersions. 
We assign the additional weight in the region k > 7i/2 to processes involving one holon 
and more than one BA spinon. In fact, that portion resembles the difference between the 
supports for one-holon/one-spinon and one-holon/three-spinon processes in the ISE model 
0. In our case, no limitation on the possible number of spinous exists, such that in principle 
all odd number of them are allowed. It is interesting to notice that using a fermionic spinon 
one is able to describe both the case J = and J = 2t. In the first case, the spinon in the 
exact solution is a fermion. At the SuSy point it is expected to be a semion 0,^] and on the 
basis of our results, we conclude that the fermionic spinon contains all possible states with 
an odd number of semionic spinous. 

Finally, we consider the quasiparticle residue Zk = {'^o^^\cj^^\'^q) at k = 7r/2 for 
J = 2t. 
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FIG. 3. Quasiparticle weight at k = tt/2 



(a) G{tt/2,-t) = G(7r/2,-r)exp 



-T (^E^ - Eq~^ (vr/2)) versus r/t. At r/t > 1 this quantity 
converges to the quasiparticle weight Z(7r/2). (b) Finite size scaling of Z{tt/2) as 

obtained from (a). The solid line is a least square fit to the form L~^^'^. We consider PJ = 30 for 
L < 48 and /3J = 60 for L > 48, to guarantee convergence in r. 
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As Fig. 2.C shows, the lower edge is very sharp and without prior knowledge, the question 
may arise whether we are dealing with a quasiparticle. is related to the imaginary time 
Green function through: 



Fig. 3. a shows G{tt/2, —t) exp (—t{Eq — Eq ^(7r/2))] versus r, where the energy difference 



is obtained by fitting the tail of G{tt/2, — r) to a single exponential form, for several sizes. The 
thus estimated Z{ti/2) is plotted versus system size in Fig. 3.b. Our results are consistent 
with a vanishing quasiparticle weight Z{tj:/2) oc L^^/^ which is the scaling obtained by a 
combination of bosonization and conformal field theory 0. Since the CPU-times scales as 
V(3 {y is the volume) the determination of the Z-factor may be efficiently extended to higher 
dimensions, in contrast to determinantal algorithms for the Hubbard model that scale as 



In summary, we have developed a new QMC algorithm which allows the determination 
of single-hole dynamics in quantum antiferromagnets. This algorithm is extremely powerful 
in the sense, that the required CPU time scales as V(3. For the one dimensional case, we 
showed that the spectral function is well described by a simple model with free spinous and 
holons with dispersions given by J and 2t respectively. The comparison of our results at the 
supersymmetric point lead to a characterization of the excitation content of the spectra for 
this particular parameter, where additional information is available from the Bethe Ansatz 
solution. Finally we computed the quasiparticle weight and showed that it vanishes as L~^/^. 
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above institutions for their support. 
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